Two-stream instability in quasi one-dimensional Bose-Einstein condensates 
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We apply a kinetic model to the predict the existence of a new instability mechanism in elongated 
Bose-Einstein condensates. Our kinetic description, based on the Wigner formalism, is employed to 
highlight the existence of unstable Bogoliubov waves that may be excited in the counter-propagation 
configuration. We identify a dimensionless parameter, the Mach number at T = 0, that tunes 
different regimes of stability. We also estimate the magnitude of the main parameters at which 
two-stream instability is expected to be observed under typical experimental conditions. 



Several aspects of BEC physics have been explored in 
recent years, with particular relevance to collective pro- 
cesses, such as collective atomic recoil lasing (CARL) 
[H, [H, the excitation of Bogoliubov waves in elongated 
condensates 0, 0, 0] , Bloch oscillations in spatially pe- 
riodic potential induced by a far-off-resonance laser field 
Q and the observation of quantum fluctuations and en- 
tanglement 

0, S, H El • Attention to the physics of BECs 
has also been paid in astrophysics, since some recent the- 
ories suggest BECs as good candidates for dark matter 
111 ]. Recent wave kinetic models have been proposed, in 
both classical and quantum regimes, to properly describe 
many cooperative phenomena in BECs. like Bogoliubov 
waves, wakefields and instabilities [IH, [H, 0, Il5| . Ex- 
act kinetic equations also seem to constitute an essential 
tool in the description of thermodynamics of condensa- 
tion. For example, Gardiner et al have applied kinetic 
methods based on the Wigner function to describe the 
formation of a BEC in the presence of a thermal bath 

In this Brief Report, we apply a quantum kinetic model 
to study the dynamical instability of Bogoliubov waves 
in two counter propagating BECs. The occurrence of in- 
stabilities in BECs has been discussed in several works, 
ranging from instabilities of the supcrflow in optical lat- 



tices [17[ and dynamical instabilities of rotating BECs 
[lit t° vortex formation mechanisms [lj)[. Recent ex- 
jerimental and theoretical works on collisions of BECs 
20l 2 1 1 . reporting on the scattering of atoms and quan- 
tum turbulence, pave the stage where such instabilities 
may take place. Furthermore, in experiments performed 
on storage ring traps for ultracold atoms [22[, where 
BECs can be split into two components which counter- 
propagate round the ring, the study of the stability crite- 
ria manifestly becomes an issue of major relevance. We 
anticipate that, in quasi one-dimensional systems, the 
splitting of the condensate into two parts that are set to 
move against each other, the Bogoliubov waves can be- 
come dynamically unstable. Such an instability should be 
referred as a "two-stream" instability, in complete anal- 
ogy with the well known phenomenon in plasma physics 
[231 ] . This dynamical instability drives an exponential 



increase of the amplitudes of the fluctuations in the con- 
densate. Such a growth induces the dephasing of the 
condensate, transferring its translational kinetic energy 
to collective and single particle excitations (phonons). 
We observe two different regimes of instability, occurring 
in both subsonic and supersonic regimes. 

In the spirit of mean-field theory, the collective field 
operator \fr(r, t) = $(r, t) + ip(r,t) of the Bose gas can 
be separated into a condensate c-number wave function 
<&(r,t) = (\I>(r, i)), where (•) stands for the expectation 
value in the ground state, and its quantum fluctuations 
ip(r,t). The condensed field verifies the normalization 
condition f $(r, t)*$(r, t)dr = N, where N is the total 
number of condensed atoms. The dynamics of the BEC 
wave function is governed by the celebrated GP equation 



iti 



~dt 



2m 



V 2 $ + (V + V SC )$, 



(1) 



where Vq = Vb(r) = m/2(w 2 _r 2 _ + uj 2 z 2 ) represents the 
asymmetric confining potential, m is the single atom 
mass, rj_ = (x 2 + y 2 ) 1 ^ 2 and z represent the transver- 
sal and longitudinal directions, respectively, and Vsc = 
<?|<i>(r, t)\ 2 is the nonlinear self-consistent potential. Here, 
g = Anh 2 a/m defines the coupling constant or interaction 
strength and a represents the s-wave scattering length. 
The main ingredient for the derivation of a quantum ki- 
netic equation is the two-point correlation function asso- 
ciated to the condensed wave function 



C(s,r) = $(r-s/2,i-T/2)$*(r + s/2,£ + T/2), (2) 

where s = r x - r 2 , r = (ri + r 2 )/2, t = (t± + t 2 )/2 and 
r = t\ — 12- By performing a double Fourier transform of 
the relative position (s) and time (r) variables, we define 
the Wigner function 



W(r, t; (J, k) = J ds J drC(s, r) cxp(icjr — ik ■ s), (3) 

where uj and k represent the frequency and the momen- 
tum of the BEC excitation, respectively. Combining the 
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GP equation |T]) with Eq^Jp]) , it is possible to derive the 
Wigncr-Moyal equation [24[ 



d h 
ih [ — — I k • V ) W = 

at m 



f da f dft r , 



(4) 



x [W- - W+] exp(-iO< + iq • r), 
where W± = W(u) ± Si/2, k± q/2) and 

V(q,n) = f dr f dtV(r,t)exp(-iq_-r + iSlt) (5) 

stands for the double Fourier transform of the total po- 
tential V = Vq + Vsc- We linearize the system by intro- 
ducing a perturbation around its equilibrium configura- 
tions, such that W± = W ± +W±, V = V Q + V, where 
W± = W±(q, O) exp(iq-r— iSlt) and V = V(q, SI) exp(iq- 
r — iSlt) represents the single Fourier components of the 
wave. We restrict our discussion to single Fourier com- 
ponents, since we are looking for the stability of a sin- 
gle mode. The same formalism can be easily extended 
to a broadband fluctuation spectrum, which is not the 
aim of the present work. Assuming that the confining 
potential is static, Vq = Vb(r), the fluctuations will be 
only rooted in the nonlinear mean-field potential. This 
statement restricts our discussion to the case of dynam- 
ical instabilities, which has no dependence on the shape 
of the external potential, in contrast to the case of en- 
ergetic instabilities f25j|. The Fourier component of the 
perturbed potential is therefore given by the convolution 
V"(q, SI) = * <6)(q, fi). Wc can easily obtain the 

linearized version of Eq. j4| 



d_ 
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^V{q, SI) [Wo+-W -]expi9, 

Z7T 



(0) 



where v = Kk/m represents the velocity field of the BEC 
and 9 = q ■ r — Sit is the wave phase. 

In order to obtain very elongated BECs, one requires 
an asymmetric trapping potential Vo(r) = Vj_(p) + V z (z), 
such that the condition = {uj 2 + u 2 ) 1 / 2 » lo z is 
satisfied. A typical experiment with 87 Rb atoms con- 
tains N ~ 10 5 atoms confined in an asymmetric har- 
monic trap with radial and axial trapping frequencies of 
~ 200 Hz and uj z /2tt ~ 20 Hz, respectively 0. 
The radial and axial Thomas-Fermi lengths of the BEC 
are R ~ 3.1 /.mi and Z m 27.1 um. Under these condi- 
tions, we can assume that the BEC is quasi-lD. Assuming 
a static density profile along the transversal direction, the 
perturbations of the self-consistent field along the z-axis 
are given by V(z,t) = V(q z , SI) eyip(iq z z — iSlt), which 
together with Eq. ^ yields 



W{z,t;w,k) = -V{z,t) — . 7 

n il — q z v z 

Integration of Eq. ([7]) over the BEC spectrum (k, co), 
together with the identity 



V(z,t) 



dk du!~ r , , . 
-^-W{z,*M (8) 



and considering that the condensate spectrum can be 
written in the form u> = u>(k), which allows the factoriza- 
tion W{z,t;w,k) = 2nW{z, t; k)5(u - w(k)), yields the 
following kinetic dispersion relation 



dk Wo(z;k-q/2)-W/ (z;k + q/2) 



(2tt) 3 



si 



q- v 



0. 

(9) 

In the one-dimensional approximation, the waves can 
not propagate in the radial direction, and we can there- 
fore factorize the equilibrium Wigncr function into its 
longitudinal and perpendicular components, such that 
W {z;k) = {2ir) 2 W {z;k z )5(kj_). Dropping the sub- 
script z, for the sake of simplicity, we can finally write 
down the dispersion relation for the longitudinal oscilla- 
tions 



9., 



qv 



Sl_ 



qv 



= 0, (10) 



where f2± = Sl±hq 2 /2m. One can easily obtain the usual 
excitation spectrum at T = 0, by simply setting the cor- 
responding equilibrium profile Wo{z;k z ) = 2imo8(k z — 
ko), which yields (SI — qv ) 2 = u 2 B q 2 + h 2 q i /Am 2 , where 
vq = hko/m is the stream speed, Ub = sj gno /m repre- 
sents the Bogoliubov sound speed and no = (n-o( r )) is 
the average background density. We immediately con- 
clude that our kinetic approach is totally equivalent to 
the usual Bogoliubov-de Gennes method. Equation IT01 
shows that these (dopplcr shifted) second-sound waves 
are dynamically stable, since SI is always real. One should 
notice that the dynamical instability has no relation with 
the Landau criterion of superfluidity (energetic instabil- 
ity), which strongly depends on the external potential. 
Mcnotti et al have explicitly shown the difference be- 
tween energetic and dynamical instab ility in the case of 
a BEC confined in an optical lattice. [25j |. 

In the counter-propagation configuration, a discus- 
sion about the stability criteria of the excitations de- 
serve some attention. In order to give a first insight 
to the problem, we assume a quasi-equilibrium profile 
describing two incoherent BEC beams, which is approx- 
imately given by Wo(r; k) = 7m [S(k — k ) + 5(k + ko)], 
where no is the mean density. This assumption is 
made by noticing that, at T = 0, the rms velocity 
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Vrms = Ti\J (k 2 ) /m of each beam is very small, such that 
QVrms <C f2, which is satisfied for modes such that l/q is 
much less than the Thomas-Fermi length Z = y/ muj z /h. 
We test the consistence of this approximation at the end 
of this work. Under these conditions, the dispersion rela- 
tion for the sound waves in the two-stream configuration 
is finally given by 



K 



-2 



1 



(3 2 (n + K) 2 -K A (3 2 (Sl~KY 



K 1 



> =0, 

( n ) 

where we have defined the normalized quantities = 
fl/tuo, K = hq/muB, and the dimcnsionlcss parameter 
(3 = vq/ub- Here, u>o = itsfco represents the 'resonance' 
frequency. The parameter (3 defines a sonic number, the 
Mach number at T = 0. It measures if the beam flow 
is either subsonic ((3 < 1) or supersonic (/3 > 1). The 
factorization of Eq. (fTTj) provides two branches 



n 2 ± 



K 

W 2 



'2 r 



1 + 2(3 2 + 2K 2 ± >/l + 8/3 2 + 16(3 2 K 2 ) 



(12) 

one of which (f2+) is always positive and describes sta- 
ble oscillations. However, the solution f2?_ is not pos- 
itive definite and negative solutions are found for the 
modes K that verify the condition 1 + 2K 2 + 2/3 2 < 
(1 + 8/3 2 + IQptK 2 ) 1 / 2 . This condition defines two dy- 
namically unstable regions in the (K, (3) plan, as illus- 
trated in Fig. ([T]). In the subsonic regime, (3 < 1, all 
the modes satisfying K < (3 are unstable; in the super- 
sonic regime, (3 > 1, unstable modes are obtained for 
(3 2 — 1 < K < (3. The representation of the dimcn- 
sionless growth rate T(K, i3) = T/ui , where T — 3(S1), 
shows that such unstable oscillations occur in both su- 
personic and subsonic regimes (see Fig.©). In the sub- 
sonic regime, all wave modes K are unstable, until the 
cut-off given by the condition K = (3. This means that 
all (?-modes are unstable up to the cut-off wave vector 
<Zc = fcoi when this mode has the same energy of the 
beam. In that case, all the wave energy is transferred 
to the condensate, and the instability vanishes. In the 
supersonic regime, (3 > 1, the picture changes. First, 
to become unstable, the phonon must have a momen- 
tum higher than the resonant value q r = "/mus/h, where 
7 = [(vq/ub) 2 — l] 1 / 2 . It means that the phonon must 
be resonant with that of the supersonic beam, so en- 
ergy can be transferred from the condensate to the wave. 
Our calculations also indicate that the increasing of (3 
increases both the value of the cut-off mode q c and the 
value of the resonant mode q r . In the subsonic (super- 
sonic) regime, the value of the maximum growth rate 
increases (decreases) as (3 increases. Therefore, the most 
unstable mode is obtained for (3 = 1 and corresponds to a 
maximum growth rate of T max ss 0.185wo, relative to the 
mode q m ax ~ 0.795gs, where qs — itiub/Ti- For the typ- 



ical experimental conditions performed with 87 Rb [261 ] . 
the average density is no ~ 10 11 cm -3 and the s-wave 
scattering length is a ~ 1 nm, as it is known for alkali 
metals 27[, which provides a Bogoliubov speed of the 
range of ub ~ 10 cm/s. Preparing the condensate beams 
at f3 = 1, we expect to observe a maximum growth pe- 
riod of t = 2T:/T max ~ 3 ms. The typical wave length 
for which two-stream instability is expected to occur is 
at the order of A = 2?t /qs ~ 1 /im, which is consistent 
with the quasi-one dimensional approximation we have 
performed. One should also notice that the validity cri- 
terion for the Bogoliubov theory \/noa 3 ~ 10~ 3 <C 1 
is also satisfied for the present experimental conditions. 
A final remark about the validity of the present results 
is in order. First, it is well known that under condi- 
tions of reduced effective dimensionality, phase fluctua- 
tions play an important role (our present discussion is 
limited to weakly interacting systems whose kinematics 
is low dimensional, but the scattering can still be seen 
as a three dimensional process, which is far away from 
the Tonks-Gireardeau regime). For such a system, be- 
low the degeneracy temperature T c i = Nhtut/ks, where 
uj t = \f^\ + w z; (the low-dimensional analogue of the 
critical temperature), one can define a second charac- 
teristic temperature, the 'phase fluctuation' temperature 
T p h = Tdhcut/fJ-, where fi is the Thomas- Fermi chemical 
potential j28|. At T ^ T p h, the condensate phase fluc- 
tuates at scales much smaller than the Thomas-Fermi 
length, which correspond to the range of temperatures 
where the two-stream instability mechanism presented 
here may be observed. Second, in the presence of a ther- 
mal gas, a more detailed calculation should include the 
dynamics of the decohcrence and the respective equilib- 
rium distribution function would correspond to that of a 
low-dimensional matter wave interferometer [29I ] . A sim- 
ple model for the condensate at finite temperature intro- 
duces a broadening of the equilibirum profile [30l ]. and 
therefore some wave modes are Landau damped, which 
suggests that our results, based on the two-independent 
beams configuration, correspond to the "best case sce- 
narium" for this dynamical instability. 

In conclusion, we have theoretically established the cri- 
teria for dynamical instability of Bogoliubov waves in 
BECs prepared in the two-stream configuration. Thresh- 
old conditions for the occurrence of unstable regimes, and 
the corresponding growth rates were derived. The condi- 
tions for the two-stream instability change in the subsonic 
and supersonic limits, where in the later the instability 
can not be excited by an arbitrarily small wave vector. 
The maximum growth rate of the dynamical instability is 
observed when the nominal velocity of the wave packets 
equals the sound speed in the condensate. The authors 
consider that the features of the instability mechanism 
presented here can eventually provide the basis for fu- 
ture experiments. 
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FIG. 1: (Color Online) Stability diagram for the two-stream 
BEC. The shadowed area represents the set of parameters for 
which the Bogoliubov waves are unstable. The upper curve 
corresponds to A' = /3. The dashed lower curve corresponds 
to the case K — \J ft 1 — 1 
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FIG. 2: (Color Online) Normalized wave growth rate T = 
Im(f!)/uo for different values of f3. Blue full line, f3 = 0.5, 
black dashed line, (3 = 1.0, and red dotted line, j3 — 1.5. The 
maximum growth rate, corresponding to the most unstable 
mode, occurs for /? = 1 and K w 0.8. 
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